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SUPPLEMENTAL DECLARATION OF PETER D. BREWER 
PURSUANT TO 37 C.F.R. § 1.131 OR. IN THE ALTERNATIVE, 37 C.F.R. 6 1.132 

Mail Stop Amendment 
Commissioner for Patents 
P.O. Box 1450 

Alexandria, VA 22313-1450a 

I, Peter D. Brewer, declare and say: 

1. I am the inventor named in the above-identified application. 

2. I refer to and incorporate by reference my earlier declaration pursuant to 
37 C.F.R. 1.131 executed on September 5, 2006 and filed in the above-identified 
application. That earlier declaration stated that I had completed the invention disclosed 
and claimed in the above-identified application in the United States of America no later 
than November 8, 2001. This supplemental declaration further explains the invention 
that was completed no later than November 8, 2001. In addition, this supplemental 
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declaration updates the declaration executed on September 5, 2006 in reference to the 
amendment to claim 1 presented in the response filed with this supplemental 
declaration. 

3. As noted in the declaration executed on September 5, 2006, Exhibit A to 
that declaration was a copy of a "Progress, Status and Management Report" for 
research on " Antimonide Based Compound Semiconductors (ABCS)" that I helped to 
prepare. The date of the report was stated on its cover: November 8, 2001. Irrelevant 
material was been redacted from the copy of the report attached as Exhibit A. The 
report evidenced actual reduction to practice of the invention claimed in this 
application before November 8, 2001. 

4. The "Description of Progress" on page two of the "Progress, Status and 
Management Report" attached as Exhibit A to the declaration executed on September 5, 
2006 contained a section that I wrote. It described research directed to a "substrate 
transfer technology focused on the preparation of the MBE grown epi-layer surfaces 
prior to wafer bonding and processes for selectively removing GaSb substrates after 
wafer bonding." "Morphological growth-defects on the surface of the Sb-based 
epilayers" are identified as problems because these defects "interfere with the bonding 
of the GaSb epilayers and the sapphire substrates." The solution was "a self-masking 
process for removing growth defects form the surface of the MBE grown Sb-based 
epilayers." 

5. As stated in the "Progress, Status and Management Report" attached as 
Exhibit A to the declaration executed on September 5, 2006, the "self-masking process" 
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involves four processing steps: "1) coating the surface of the wafer with a thick 
photoresist layer (5-10 microns), 2) dry-etching the resist layer to a thickness of -0.5 
microns (to reveal the tops of the defect structures but protecting the remainder of the 
semiconductor surface), 3) wet chemical etching of the exposed defect structures, and 4) 
stripping of the remaining photoresist layer." 

6. The "Progress, Status and Management Report" attached as Exhibit A 
observed that "this process effectively removes the protruding defect structures from 
the surface of the semiconductor wafer without [affecting] the surrounding epilayer 
material." The "Report" noted the successful results of the process: "[i]nitial results 
using this process to prepare as-grown HBT wafers for bonding to sapphire substrates 
indicate enhanced bonding yields as a result of eliminating the morphological growth 
defects. In these experiments, bonding surface area yields as high as 94% were 
obtained." 

7. The step of "1) coating the surface of the wafer with a thick photoresist 
layer (5-10 microns)" will inherently result in the formation of a layer of photoresist on 
the surface of the wafer that will have a planar surface on its top surface, that is, the 
surface not in contact with the wafer. The reason for this is the thick photoresist will be 
planarized by its own surface tension and viscosity. The defects will not affect the 
planarity of the surface of the photoresist layer because the layer is much thicker than 
the defects. Note that the second step is "2) dry-etching the resist layer to a thickness of 
-0.5 microns (to reveal the tops of the defect structures but protecting the remainder of 
the semiconductor surface)," which indicates that the defect structures will have heights 
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above the general surface of the wafer of no more than about one-tenth to one twentieth 
of the thickness of the thick photoresist layer. 

8. I attach as Exhibits B and C to this supplemental declaration copies, 
respectively, of two articles: Jun-Bo Yoon, et al, Planarization and Trench Filing on Sever 
Surface Topography with Thick Photoresist for MEMS, SPIE Vol. 3511 (September 1998), at 
pp. 297-306 ("Yoon, et al/') and Peter C. Sukanek, A Model for Spin Coating with 
Topography, 136 (No. 10) J. Electrochemical Soc. 3019-3026 (October 1989)("Sukanek"). 

9. Yoon, et al. (Exhibit B) shows that to achieve a value of planarity (/3) 
greater than 95% (0.95) the value of the thickness of the photoresist (do) relative to the 
topography of the surface (to) the value of do/to must be greater than about six (6) for 
AZ4562 photoresist. See Figure 6(b) of Yoon, et al. at p. 303. The defect structures 
discussed in my Exhibit A have heights of between one and ten microns. Persons of 
skill in the art will understand that typical defects are crystallographic defects (also 
known as oval defects) and will have a height of about one micron and that defects with 
a greater height are only occasionally observed and are easily located and removed by 
known processes. Six times the height of the typical defect is 6.0 microns. A "thick 
photoresist layer (5-10 microns)" therefore will have good planarity. The planarity will 
be greater than 90%. 

10. Sukanek (Exhibit B) discusses some of the physical factors that result in 
greater planarity as a result of using thick photoresist films. Sukanek assumes (at page 
3020, item h) that the characteristic height of the substrate variations Op (that is, the 
topography of the defects on the surface) is much less than L (the thickness of the 
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photoresist). Sukanek states in his abstract that the "viscosity and surface tension of the 
film [of photoresist], together with their variation with solvent content, play an 
important role in determining whether planar or conformal surfaces are achieved. 
Surface tension tends to planarize the film, whereas surface tension gradients tend to 
make the film more conformal. If viscosity rises very quickly as a result of solvent 
evaporation, the film topography becomes 'frozen-in, 7 and surface forces cannot 
planarize/ 7 This behavior is described mathematically in equation 39 on page 3022 and 
is described in words in the paragraph below it: "[f]rom this result, it s seen that a 
planar surface (yn S =0) is achieved . . . when the surface tension is very large/ 7 

11. Accordingly, a person of skill would know that the step of y/ l) coating the 
surface of the wafer with a thick photoresist layer (5-10 microns) 77 would inherently 
result in the formation of a layer of photoresist on the surface of the wafer that will have 
a planar surface on the top surface of the photoresist 

12. I understand that claim 1 of this application is being amended to read as 
follows: 

A method for removing defects from a semiconductor surface, 
comprising: 

coating the semiconductor surface and the defects with a protective 
layer, wherein the protective layer has a planar top surface; 

thinning the protective layer to selectively reveal top portions of 
the defects; 

removing the defects; and 

removing the protective layer. 
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This method is disclosed in Exhibit A, the "Progress, Status and Management Report/' 
The "protective layer" reads on the "thick photoresist layer," which as mentioned 
above, inherently had a planar top surface because of its own surface tension and 
viscosity. 

I declare further that all statements made herein of my Own knowledge are true; 
that all statements made herein on information and belief are believed to be true; and 
further that these statements were made with the knowledge that willful false 
statements and the like are punishable by fine or imprisonment, or both, under Section 
1001 of Title 18 of the United States Code and that such willful false statements may 
jeopardize the validity of this application or any patents issuing thereon. 



Date; May 3, 2007 




Peter D. Brewer 
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Description of progress 

Work included establishing a wafer transfer process, and improvement of the 
performance of existing prototype ABCS-based HEMT and HBT devices. 
Substrate Pevelonmentt 

The development of tbe substrate transfer technology focused on the preparation 
of the MBE grown epi-layer surfaces prior to wafer bonding and processes for selectively 
removing OaSb substrates after wafer bonding- 

As reported in the previous progress report, morphological growth-defects on the 
surface of the Sb-based epilayers are found to interfere with the bonding oi the GaSb 
epilayers and the sapphire substrates. The growth-defects and other partioleH found on 
the wafer surfaces, cause the bonding wafers to deform around them forming, circularly 
un-bonded interface areas or voids. The un-bonded areas resulting from oven small 
protuberances are fairly large, for example, a particle of -1 micron diameter leads to an 
un-bonded area with a diameter of -0.5 era for GaSb and sapphire wafers. The defects 
are produced during the MBE growth of epilayers and are the result of el fusion cell 
spitting or growth defects caused by impurities or surface imperfections on the original 
GaSb substrate wafer. The defects are an integral part of the semiconductor wafer 
surface and cannot be removed with conventional particulate removal processes. The 
density of defects on the wafer range from M00/cm 2 and range in size from 1-50 
microns with heights typically from 1-10 microns. 

HRL has developed a self-masking process for removing growth-defects from the 
surface of the MBE grown Sb-based epilayers. This process involves four processing 
steps: 1) coating the surface of the wafer with a thick photoresist layer (5-10 microns), 2) 
dry etching the resist layer to a thickness of -0.5 microns (to reveal the tops of the defect 
structures but protecting the remainder of the semiconductor surface), 3) wet chemical 
etching of the exposed defect structures, and 4) stripping of the remaining photoresist 
layer. Although simple, this process effectively removes the protruding defect structures 
from the surface of the semiconductor wafer without effect title surrounding epilayer 
material. Initial results using this process to prepare as-grown HBT wafers for bonding 
to sapphire substrates indicate enhanced bonding yields as a result of diminating the 
morphological growth-defects. In these experiments, bonding surface area yields as high 
as 94% were obtained. 

A second area of this task focused on developing a process for selectively 
removing GaSb substrates from InAs epilayers after wafer bonding, Tlus month's 
activity has centered on developing a three-step process to thin and selectively remove 
the GaSb substrate. The process includes the following steps: 1) lap and polish the GaSb 
substrate (bonded to sapphire) to a thickness of 50 microns, 2) continue thinning of the 
GaSb substrate material to ~20 microns using a selective wet chemical etch, and 3) 
remove thinned GaSb material from the InAs epilayers using a highly selective dry etch 
process. The multi-step process is used to efficiently remove the substrate material from 
the epilayer device structure and not over burden the dry etch system with excessively 
long etch runs> The dry etch process uses an Inductively coupled plasma etch system 
using a Cl-based etch chemistry mat has demonstrated etch selectively oi' 1000:1 for 
GaSb and InAs, respectively. This process was successful used to completely remove the 
GaSb substrate material from a InAs diode structure epilayer bonded to sapphire. 
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Planarization and trench filling on severe surface topography with thick 

photoresist for MEMS 



Jun-Bo Yoon, Gilbert Y. Oh, Chul-Hi Han, Euisik Yoon, and Choong-Ki Kim 

Department of Electrical Engineering, KAIST, 
373-1 Kusong-dong, Yusong-gu, Taejon 305-701, Korea 

ABSTRACT 

We have examined a simple and low-cost method to achieve planarization and trench filling on a severe surface topography 
for MEMS. The method simply uses a single-layer coating of a thick photoresist or polyimide, where the coating thickness 
is much greater than the severe surface topography. 

From extensive experiments, we extracted simple empirical formulae for the planarization factor (3 of die thick photoresist 
AZ4562 and polyimide PI261 1, which let us know the minimum film thickness required to obtain a certain P on a given 
surface step height and pattern density. We could compare the planarization capability of AZ4562 and PI2611 by this 
method. 

Moreover, after the planarization with a thick photoresist, we have shown a new way, other than the plasma etching, to 
remove the upper photoresist layer conformally. Therefore, we could remain the photoresist only in trenches and fill up very 
deep and wide trenches with the photoresist. Also, we have shown the etch-back result of PI261 1 using conventional 0 2 
plasma REE. 

Using these methods, we obtained the P of 98% on the surface of 20um-deep and 200um-wide lines and spaces with a 
single-coated 70um-thick photoresist And lOum-deep and 200um-wide trenches as well as 2um-deep and 50um-wide 
trenches were neatiy filled up with the photoresist As applications for MEMS, We fabricated microchannels by metal 
coating after the trench fill-up process, and even multilevel microchannels by controlling the height of the photoresist 
remaining in the trench and repeating the trench fill-up process. 

Keywords: Planarization, etch back, trench fill, thick photoresist, polyimide, microchannel, MEMS, planarization factor 

1. INTRODUCTION 

Planarization is the term that becomes more and more important along with advances in IC technology in the ULSI era and 
emerging MEMS. As an example, multilevel metallization introduces recent Chemical Mechanical Polishing (CMP) for 
global planarization in VLSI application. 1 For MEMS, the topography of the substrate is increasingly modified: bulk 
rmcrormchining etches the substrate, surface rmcrormcriining usually deposit sacrificial layer and high-aspect-ratio 
microstructures on the surface. Once the surface of the substrate has severe nonplanarity, usually no more process can be 
possible due to the limited depth of focus of the exposure optics in projection aligner and gap between mask and wafer in 
contact aligner. 

There are many techniques for planarizing surface with topography, mainly three methods: sacrificial layer and etch-back 
(EB), 2 in other words, plasma planarization, 3 selective dielectric deposition (SDD), 4 and CMP. 1 In SDD, surface condition 
can easily affect the dielectric deposition, and for CMP, global planarization can be possible, however, the method suffers 
from a dishing effect, surface scratch, and residual contamination. 
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In EB, a low viscosity sacrificial layer (usually photoresist or polyimide) is used to spin-coat the uneven substrate surface to 
obtain a flat surface. Subsequently, plasma etching is used to etch-back the sacrificial layer with equal etch rate with the 
substrate as shown in Figure 1. Because it is simple and cost-effective, many types of approaches have been studying for 
EB: two or multi-layer planarization (TLP or MLP), 2 * " image reversal and coating (IRC), 8 * 9 defocused resist patterning with 
blanket stripe mask (DRESS), 10 and planarization block mask (PBM). 11 These need additional lithographic steps or 
photomasks, or suffer from interactions between multiple polymers. By the way, planarization with single polymer layer 
was attempted, 12 however, the thickness of the sacrificial polymer layer was limited about one micrometer in that time, so 
that the degree of the planarization was not good. 

In this paper, our key idea is very simple that applying the film much thicker than the surface topography. Therefore, we 
have approached the single-layer planarization with thick photoresist and polyimide. Using single or multiple coatings of 
the same photoresist with thickness over lOO^im, if necessary, we have shown successful planarization results on the surface 
with severe topography for MEMS as well as IC applications. And from extensive experiments, we have extracted a simple 
empirical model, which let us know the minimum film thickness required to obtain a certain planarization level on a given 
surface step and pattern density. 

Moreover, after the planarization, we have experimented an easy way to remove conformally the upper photoresist layer to 
remain the photoresist only in trenches, and hence to fill up very deep and wide trenches with the photoresist. Also, a 
plasma etch-back experiment for polyimide has been shown. 

As an application of this technique, we have fabricated microchannels by metal coating after the trench fill-up process, and 
even multilevel microchannels by controlling the height of the photoresist remaining in the trench and repeating the trench 
fill-up process. 

2. PLANARIZATTONWITH THICK PHOTORESIST 

2.1 THICK PHOTORESIST AZ4562 AND POLYIMIDE PI2611 

Thick photoresist AZ4562, the thickest one in AZ4000 series manufactured by Hoechst, is a Novolac type, highly viscous 
(400cSt), and highly transparent photoresist It has been frequendy used for MEMS applications: TAB/bump process, 13 
high-aspect-ratio photolithography, 14 and 3D UV-microfonning. 15 Mainly, this photoresist is used as a high-aspect-ratio 
mold for electroplating. 

To obtain the relationship between the film thickness and the spin speed as well as the spin time, we coated the photoresist 
on 4-inch wafers by a conventional spin coater and measured the film thickness. Figure 2 shows the result. Figure 2(a), 
which was presented by the vendor, shows the film thickness, dO, is proportional to the inverse square root of the spin rpm, 
and Figure 2(b), which was own experimented, also shows the film thickness is proportional to the inverse square root of die 
spin time. 

We obtained the film thickness over lO^m by reducing the spin time, not by reducing the spin speed, as the vendor 
recommended for the film uniformity. For the film thickness uniformity, l^m standard deviation in 50\im film thickness 16 
and uniformity variation with the spin speed was reported elsewhere. 13 In our experiment, uniformity of ±5% for single 
coating of the 90|im-thick photoresist over a 4-inch wafer was observed excluding the thick edge bead. For thicker films, 
multiple coating is effective. No interaction was observed between layers with a sufficient soft baking, which is consistent 
with the fact of double thickness with double coatings. 9 

For the supplementary and comparison purpose, polyimide PI261 1, manufactured by DuPont, has also been used. Figure 3 
shows the relationship between the film thickness and die spin speed as well as the spin time of PI261 1 . 

2.2 EXPERIMENTS 

2.2.1 DEFINITION OF THE PLANARIZATION FACTOR p 
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Prior to describing the experiment procedure, the term planarization factor p needs to be defined as it pertains to our work. 
As shown in Figure 4(a), the definition of the planarization factor p is given by, 17 



^ttep 



where t' step and t r step are the initial and final step heights, respectively. In complete planarization cases, p=l; if no 
planarization exists, then P=0. And according to Figure 4(b), our |3 can be expressed as, 

_ , (d2+t0)-d\ 2dl-d2-t0 

P=l 7l = ~ . (2) 

where tO, dl, and d2 are the surface step height, film thickness in the groove, and film thickness on the mesa, respectively. 
2.2.2 TEST VEHICLE PREPARATION 

For making surface steps with vertical cross-section, we used anisotropic etching of (110) surface-oriented silicon wafer. 
The surface step height varied from 2^m up to 20nm, and pattern density from 5^im lines & spaces (L&S) up to 200^im 
L&S. The L&S pattern had SOOO^im in length. The KOH-etched wafer was RCA-cleaned and HMDS vapor-primed for 
photoresist coating. We fixed the spin speed to 2000 rpm and varied the spin time to obtain the film thickness from 10^m up 
to lOO^im. For the film thickness over 40^m, we used multiple coating. After all the single coating of the photoresist, the air 
stabilization for 20min is followed by the soft baking at 90°C for 20min on a hot plate. 

For polyimide PI261 1, L&S patterns were fabricated on the flat wafer by electroplating copper up to 10 jam thick with the 
same photomask used for the silicon etching. Over the copper L&S pattern, PI2611 was coated at 2000ipm with various 
times, then soft-cured in convection oven with N 2 ambient at 200°C for 30min followed by hard curing at 300°C for 30min 
in the same environment 

2.23 THICKNESS MEASUREMENT 

Samples are cleaved and all thickness parameters are obtained from viewing the cross-section by the ruler-armed optical 
microscope. The cleaved plane is (111) plane, and is perpendicular to the (110) wafer surface. For PI2611, the surface 
profile was measured using a-step manufactured by Tencor. 

23 PLANARIZATION RESULTS 

Figure 5(a) shows the 20nm-depth step with 200nm L&S, which was planarized by 70nm-thick photoresist to have a p of 
98%. And Figure 5(b) is for the 20nm-depth step with 20^m lines and 90nm spaces, which was planarized by 23^im-thick 
photoresist to have a p of 93%. The Figure 5(c) was obtained after applying 30nm-thick PI2611 on 9.5^m-thick copper 
patterns of lOOjim L&S, and t f stcp of only lfim was observed (P=96%). 

2.4 EMPERICAL FOMULA FOR p 

We extensively experimented for planarization by applying films of various thickness on various surface topography. The 
planarization factor P was surely affected by the step height, pattern density, and film thickness. There already has been an 
empirical model for calculating dl, the photoresist thickness over space as given by, 18 

dl = dO + |° jf * /r(r,G )f(r/s)rdQdr 9 (4) 

where h(r,9) is the actual height of the surface topography; r and 9 are the cylindrical coordinates, f is a monotonically 
decreasing weighting function of r and s, e.g., exp(-r/s), with s being the characteristic length of the material. 

However, we needed a more easy and practical formula to know the mmi'mnni film thickness required to obtain a certain p 
on a given step height and pattern density. We tried to extract an empirical formula for p using an experimental result that 
when dO and d2 were normalized to the surface step height (dO/tO, d2/t0), many experiments showed linear relationships 
between dO/tO and d2/t0. 19 They just used the result for observing photoresist step coverage, however, we used it for 
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estimating the |J. And we observed that the normalized dl also showed linear relationships with dO/tO as shown in Figure 
6(a), where die data were extracted from the 2(%m L&S pattern with 2~20jim step and KMOOjim film thickness. From 
the linear fitting parameters, die formula for (J can be expressed as, 





P = 




(2al-a2)— + 2M-62-1 




(3) 



dl 
tO 



y 



where al, bl, a2, and b2 are the linear fitting parameters. Therefore, from several test measurements, we can extract the 
parameters for a given pattern density, and then, we can have an empirical formula for P with a function of only dO/tO. 

Figure 6 shows the linear fitting parameters, empirical formulae, and measured P of A24562 for the L&S of 200|im/200|im. 
And Figure 7 shows the similar results for PI261 1. From Figure 6 and 7, we can see the linearity is well observed. In Figure 
8, we can see the relationship between p and pattern density, and we can compare the planarization capability of AZ4562 



From these empirical formulae for p, we can estimate how much thick the applying film should be to planarize lO^m-thick 
300|im lines and 100 jam spaces with the P of 80%. The answer can be easily obtained from Figure 8 that a 20 |im- thick film 
of AZ4562 or 40jim-thick film of PI261 1 will be more than enough, because the 200^m L&S data covers the case of 300^un 
lines and 1 OOfim spaces. 



In addition to planarization, we investigated etch-back methods. For thick photoresist AZ4562, we tried a new method, not 
the conventional gas-phase plasma etching, that is Underexposure and Overdevelopment Technique (UOT). Figure 9 
illustrates this simple concept Once the surface is planarized by the thick photoresist, the photoresist is to be exposed with 
no photomask and insufficient exposure dose. Then the development is to be proceeded until only the trench is filled up 
with the photoresist. Figure 10 shows the development rate is notably decreased as the development proceeds, when the 
exposure dose is not sufficient, therefore, the remaining thickness is well controlled by the exposure dose. Also, UOT can 
be repeated several times to etch very thick film, and used to pre-etch the thick film to a certain thickness level that is 
affordable to the following plasma etching step. 

We obtained several results on the trench filling as shown in Figure 11. The V-grooves in Figure 11(a) has the depth of 
15|im, the trench in Figure 1 1(b) has the width and depth of 200nm and lOum, respectively, and the trench in Figure 1 1(c) 
has the width and depth of 50jim and 2.3nm, respectively. 

After the trench filling, metal layer can be vacuum-deposited (and electroplated, if thicker layer is preferable), then the 
photoresist is removed to make microchannels for MEMS. Figure 12(a) shows 40^m-wide, 2.3^m-deep microchannels. 
And Figure 12(b) shows lOjim-wide, 2.3nm-deep, and 5jim-separated microchannels that are corresponding to 1693dpi 
nozzles for an inkjet printhead. Also, Figure 12(c) shows 6jim-wide, 7.8nm-deep, and 4^m-separated microchannels that 
are corresponding to 2540dpi nozzles. Also we fabricated multilevel microchannels by repeating UOT and metal deposition 
along with controlling photoresist depth in the trench during UOT, and the fabricated two-lever microchannels are given in 
Figure 13. 

In addition to UOT, we applied the conventional plasma etching to etch-back of PI261 1. We used the Reactive Ion Etching 
(RIE) with the condition of 100W, 150mTorr, and 20sccm 0 2 . The etching rate of PI261 1 was 0.34nm/min. We etched the 
sample of Figure 6, which was 30^m-thick PI261 1 on 9.5jim-thick copper patterns of lOO^rn L&S, and obtained the surface 
profile as shown in Figure 14(a) after 40min etching. We got Figure 14(b) just after die etch front met the copper, we 
stopped etching when we observed the appearance of copper by the viewing port in the RIE. Figure 14 together with Figure 
6 shows that the planarized surface profile was maintained during die RIE. 
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4. SUMMARY 



In this work, we applied thick photoresist and polyimide to planarization and trench filling on a severe surface topography 
for MEMS. The l(M00um-thick film was tried to planarize surface steps of up to 20um-deep and 200um-wide L&S. 

From the empirical formula of the planarization factor, we could estimate the minimum film thickness that is required to 
obtain a certain planarization factor on a given surface step and pattern density. And we could compare the planarization 
capability of the photoresist AZ4562 and the polyimide PI26 1 1 by this formula. 

Also, we could fill severe trenches with the photoresist by a new simple method, UOT after the planarization- For etch-back 
of the polyimide, we used the conventional plasma etching with 0 2 plasma RIE and observed good preservation of the 
surface profile. 

These methods are very simple, therefore cost-effective, and adequate to planarize and fill the severe surface topography 
usually occurred in MEMS. As applications for MEMS, We fabricated microchannels and even multilevel microchannels by 
metal coating after the trench fill-up process followed by removing the photoresist. We showed microchannels 
corresponding to 2540dpi nozzles for an inkjet printhead. This technique is expected to be easily applied to other 
microfluidic device fabrication as well as ultra-high density inkjet printheads. 
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Figure 1. Planarization by sacrificial layer and etch-back. 
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Figure 2. The film thickness of AZ4562 with coating conditions, (a) Spin speed vs. 
film thickness (b) Spin time vs. film thickness 



15 



£10 

c 

o 
'£ 

~ 5 



Spin time : 30s 



After prebake 140°G, 
3min, hotplate 

/ 




26966 



rpm 



1 684 1 



rpm 

After cure 350°C 30min, H 



2000 



3000 4000 
Spin speed (rpm) 



5000 



40 



30 



0) 

c 



20 



CD 

rT 10 



Spin speed : 2000rpm 



After soft cure 200°C, 
\y Am 30min, N, 




A filer hard cure 
300;C, 30min, N 2 



d0 = 



10 15 20 

Spin time (s) 



Figured, the film thickness of PI2611 with coating conditions, (a) Spin speed vs. 
film thickness (b) Spin time vs. film thickness 



(a) 



t'step 



i 



t'step 



(b) 



1111m a 

MmSmmM Planarization 



material 




p =1— ^ 



0 _ x (dl + tQ)-d\ 2d\-d2-tQ 
d\ d\ 



Figure 4. The definition of the planarization factor p. (a) From the reference 17 (b) 
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Figure 6. An empirical formula for p from the linear relationships between dl/tO and dO/tO, 
d2/t0 and d0/t(h (a) Linear parameters extraction (b) Obtained empirical formula 
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Figure 8. Comparison of the planarization factor p by empirical formulae. 
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Figure 9. Underexposure and Overdevelopment Technique (UOT). 
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Figure 10. The overdevelopment characteristic with varying the exposure time. 



00 




Groove depth : \5yim 



Figure 11. Several results on the trench filling with thick photoresist by UOT. (a) 
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Figure 11. (continued) (b) the deep and wide trench (c) the shallow and wide trench 
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Figure 12. MicroChannel made by UOT and metal deposition, (a) Cross-sectional view of the 
shallow and wide microchannels (b) Corresponding to 1693dpi (dots per inch) nozzles for an 
inkjet print head (c) Corresponding to 2540dpi nozzles 
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Figure 13. Multilevel microchannels made by repeating UOT and metal deposition, (a) Cross- 
sectional view of the two-level microchannels before removing the photoresist (b) Magnified view 
of them after removing the photoresist 
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A Model for Spin Coating with Topography 

Peter C Sukanek*' 1 

Philips Research Laboratories, 5600 J A Eindhoven, The Netherlands 

ABSTRACT 

A theoretical model for spin coating on a substrate with topography is presented. The equations assume a thin film 
wUh uniform properties in the direction normal to the substrate. By assuming a small amplitude for the substrate varia- 
tions two sets of coupled equations are derived, one for the mean value of the thickness and its composition, and another 
for the fluctuations. The viscosity and surface tension of the film, together with their variation with solvent content olav 
an important role in determining whether planar or conformal surfaces are achieved. Surface tension tends to planarize 
the film, whereas surface tension gradients tend to make the film more conformal. If viscosity rises very quickly as a result 
of solvent evaporation, the film topography becomes "frozen-in," and surface forces cannot planarize 



Spin coating is widely used in the microelectronics in- 
dustry to produce thin, uniform coatings of materials from 
solution on substrates. The most common type of material 
is an organic photoresist, which is then used in a litho- 
graphic operation to delineate patterns. In many applica- 
tions, the substrate on which the resist is spun is not uni- 
form, but processes some variation in height as a result of 
previous processing. The relation between the topography 
of the substrate and that of the resulting film is important. 

In most applications, it is desirable to achieve a planar 
resist surface, regardless of the extent of underlying height 
variations. This is the case, for example, in lithography 
with a small depth of focus, where slight variations in the 
distance of the resist film from the optical system would 
adversely affect the quality of the latent image. A planar 
film is also required in etch-back metallization. In this pro- 
cess, a metal layer is formed with "pillars," and then 
covered with a conformal dielectric. A planar resist film is 
then applied, and etched in such a way that the resist and 
dielectric are removed at the same rate. At the end of the 
process, all that remains .are the metal pillars protruding 
from the dielectric. These then form the contacts for the 
second layer metal. 

In other applications, however, it may be advantageous 
to achieve a conformal film, This is the case for the dielec- 
tric discussed above, which could be spun onto the sub- 
strate as a polyimide or spin-on glass. In lithography with a 
large depth of focus, conformal coverage would not com- 
promise the image quality and may reduce the necessity of 
over-developing the pattern in some areas. 

In a series of publications, White (1-5) has described an 
empirical model for coating on topography. The substrate 
surface is expressed as a Fourier series. The coating sur- 
face is assumed to be of the same form, but with different 
Fourier coefficients. These coefficients are determined by 
applying a low pass frequency filter to the substrate coef- 
ficients. The filter is assumed to be characterized by two 
empirical parameters which must be determined for each 
different photoresist and for different operating condi- 
tions, such as the spin speed and post-apply bake tem- 
perature. 

A different, but related, approach was taken by Wilson 
and Piacente (6, 7). These authors define an "average struc- 
ture height" in terms of an undefined function of position 
on the surface. This function is different for different types 
of resists and different processing methods. While the au- 
thors do not do so, the function can apparently be deter- 
mined empirically by measuring the resist height on wa- 
fers with different types of topography. 

Stillwagon et al (8) take a more fundamental view of the 
leveling process. These authors derive an expression for 
the "leveling time," the time required for a non-evaporat- 
ing film to planarize as a result of a surface tension-driven 
flow. The film is assumed to be conformal after the coating 
process, and then flow during a leveling period. In their ex- 
periments, liquid epoxies were used, which were then 
cured, and the resulting topography measured. Their 
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"model" was meant to indicate the relative effect of impor- 
tant variables such as mean thickness, viscosity, surface 
tension, and feature size. 

While the empirical approaches mentioned above ap- 
pear to be successful, it would be advantageous to have a 
model available which relates the film topography to the 
physical properties of the resist as well as the coating and 
post-apply conditions. For use with photoresists, solvent 
evaporation effects should be included, as well as the vari- 
ation of surface tension and viscosity with the solvent con- 
tent. Previously, a model was developed for spin coating 
on flat substrates (9). The model predictions are in good 
agreement, at least qualitatively, with available data. It 
predicts, for example, that the film height increases with 
the initial solids fraction, the initial viscosity of the solu- 
tion and with the solvent volatility, and that it is inversely 
proportional to the square root of the spin speed. In this 
paper, the model is extended to predict the variation of 
thickness on a substrate with small, periodic height varia- 
tions. 

The model draws on earlier work by Overdiep on the lev- 
eling of paints (10-1 1). In many ways, the problems are sim- 
ilar. Overdiep shows that the coating surface over a rough 
substrate depends on the surface tension of the coating 
material and how the tension varies with composition dur- 
ing the drying process. In Overdiep's model, which as- 
sumes a constant evaporation rate during the drying, ei- 
ther a planar or a conformal coating can be achieved 
depending upon whether surface tension or surface ten- 
sion gradient effects dominate the process. The pre- 
dictions of the present model agree qualitatively with 
those results. However, there are differences. In particular, 
substrate rotation must be accounted for in the present 
analysis. In addition, a decreasing evaporation rate is in- 
cluded in the present model, as is a more complete de- 
scription of the normal stress condition at the film surface. 

Theoretical 

A complete solution to the problem would be quite com- 
plicated. In the present treatment, a number of assump- 
tions are employed. The major ones are as follows: 

a) The time scale for the thinning, t\ is much larger than 
that for the rotational flow, H. Hence, for the equations 
used here, the rotational flow appears only as the driving 
force for the radial velocity. 

b) The characteristic length for the film thickness, L, is 
much smaller than the substrate radius, R. Consequently, 
the lubrication approximation may be employed. 

c) The fluid has two components, a volatile solvent and a 
nonvolatile solute. 

d) The process is isothermal. 

e) The film is uniform throughout its thickness. Hence, 
it may be characterized by a single value of the solvent 
fraction, x, which may only depend on time and the radial 
and azimuthal coordinates, r and 0. 

f) Solvent evaporation may be modeled as mass transfer 
from a rotating disk. The mass transfer coefficient is con- 
stant over the surface, and increases as the square root of 
the spin speed. 
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* g) The viscosity of the fluid depends on the liquid frac- 
tion, but is independent of the shear rate. 

h) The characteristic height of the substrate variations, 
a p , is much less than L. 

i) The substrate is periodic in r and 6. 

A critique of assumptions a) through g) is provided by 
Bornside et. al (12). The most restrictive assumption is 
probably e). As solvent diffuses from the film, the diffu- 
sivity of the material drops, making further evaporation 
more difficult This would cause axial gradients in the film, 
which are ignored here. 

Previously (9), the complete equations of conservation of 
mass and momentum were put into dimensionless form 
and simplified with assumption b). While a similar ap- 
proach could be employed here, it is simpler to make the 
standard lubrication approximations (13) at the outset. In 
addition to the radial flow driven by the rotation of the 
substrate, surface tension effects at the film-air interface 
drive an azimuthal flow, as well as affecting the radial flow. 

The simplified radial and azimuthal components of the 
equation of conservation of momentum are (symbols are 
defined in the List of Symbols section) 



dp d 2 V z 

-prH 2 = + fx — — 

dr dz 2 

1 dp a 2 i> 0 
0=- - + jt— ^- 

r ae sz z 



[1] 



[2] 



The fluid is assumed to adhere to the substrate at the sur- 
face, z D . At the top of the film, z = h t the stresses are bal- 
anced by surface tension variations 

[3] 



[4] 



These equations are easily solved for the velocity pro- 
files and the flow rates (per unit circumference) in the ra- 
dial and azimuthal directions. These are given by 
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Total and species mass balances on a differential area of 
the film give the following equations for film thickness and 
solvent fraction 



ah la la 

— = - — — (rq r )- — q B -e/p 

at r dr r ae 



[9] 



dx dx 1 ax 
(h - z D ) = -q r q 0 e(l - x)/p [10] 

at ar r ae 

where e is the solvent flux from the film resulting from 
evaporation. With assumption f), this may be written as 



and k is the mass transfer coefficient, with units of veloc- 
ity. For laminar mass transfer on a rotating disk (9) 



k = oH 0 - 5 



[12] 



Included in Eq, [1 1] is a parameter, x R , to account for sol- 
vent which remains in the film after processing. It is 
known that there is always some solvent remaining after 
spinning, and also after a typical post-apply bake (14). The 
probable physical reason for this solvent retention is the 
greatly reduced diflfusivity through the solvent-depleted 
surface regions of the film, as discussed above. The film 
does not come to equilibrium with the atmosphere under 
normal processing conditions. Since the present model as- 
sumes uniform properties throughout the film [assump- 
tion e)], this effect cannot be accounted for. Thus the resid- 
ual solvent must be considered an empirical quantity. 

The above equations must be solved with an expression 
for the pressure. This is given by a balance between pres- 
sure, viscous and surface tension forces at the film-air in- 
terface 

dv 2 f 1 a / dh \ l d 2 h 1 

-p + 2jx « a r + [13] 

dz lr dr V dr) r 2 ae 2 J 1 J 

It is convenient to put the defining equations into di- 
mensionless form using the following variables 



y = h/L ^ « i = z/L P = pL/a° 

t = tlV 7 = a/a° i = r/L u = v z t'/L 



[14] 



where L and t' are the characteristic length and time scales 
used previously (9) 



L = (ajjL°/p) ,/3 n- ,/2 

v = (kl°/ P ) ,/3 a- m a- 1 



[15] 



The ratio LJt' is the mass transfer coefficient, a characteris- 
tic velocity for the problem. In the absence of rotation 
(during baking, for example), the dimensionless distance 
and times scales are then 



L B = /i B ° 
V B = h B °/k 



[16] 



where the subscript B is used to indicate baking condi- 
tions, and o indicates the value prior to baking. 

The dimensionless equations for the film height, solvent 
fraction, and pressure are 
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•where T is the variation of surface tension with solvent 
fraction 



dX 



9" is a differential operator 



dv\ 1 d^y 



[20] 



[21] 



and Ca is the capillary number, the ratio of viscous to sur- 
face tension forces 



Ca = \L°k/a° 



[22] 



In these equations, a "rotation parameter," % has been 
introduced. Those terms which & multiplies arise because 
of the rotation of the substrate; the other terms arise from 
evaporation and surface tension effects. Hence, when the 
substrate rotates, 0t = 1. In a typical spin coating opera- 
tion, the film is baked after the apply step. During the 
bake, SftisO. 

The surface roughness is assumed to be of the form 



£ 0 = a(p(?, 6) 



[23] 



In order to solve these equations, the amplitude of the sur- 
face roughness is assumed to be small. The film height y, 
solvent concentration x, and pressure P are expanded in a 
series in terms of a 



V = Vo + *Vi + a 2 y 2 + . . . 
x = x 0 + axj + a 2 x 2 + . . . 
P = P 0 + aP x + a 2 P 2 + . . . 



[24] 



The residual solvent content, x R , is assumed to be con- 
stant. 

These expressions are substituted into Eq. [17]-[19], and 
terms of like power in a equated. At zero order (a smooth 
substrate), the height is assumed to be independent of po- 
sition. This is in agreement with previous results for spin- 
ning on a flat substrate (9). The equations for the first two 
orders of x and y are 
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The equation for the pressure at first order is 



-P«-2Ca^-^[aVi-9>] 

+ y 0 2 &Pi - 2yJT9x l = -7%, [29] 

Equations [25] and [26] are identical to those found previ- 
ously (9). The parameter N is the variation of viscosity with 
solvent fraction 



AT = 



ax 



[30] 




Fig. 1. Surface geometry indicating wavelength (X) and amplitude 
(a) of the surface roughness, the mean film height (yj, and the surface 
amplitude (y\). With y { = 1, the surface is conformal; with yi = 0, the 
surface is planar. 

Before further progress can be made in solving the equa- 
tions, the surface roughness must be specified. In what fol- 
lows, we assume that the surface has a periodic roughness 
in one direction only. In this case, the most general expres- 
sion for (p is 



9 = ly sn sin (nA£ cos 6) 



[31] 



with <p 8 i = 1. The geometry is illustrated in Fig. 1 for the 
case of sinusoidal variations. More complicated forms of 
the roughness can be easily treated by including cosine 
terms in the expansion. However, non-periodic structures 
and isolated features cannot be treated in this manner. A 
Fourier integral technique must be used in that case. 

With this simple form, the film thickness, solvent frac- 
tion, and pressure may be assumed to depend on a single 
parameter, x 



x = € cos e 



[32] 



The solution of these equations is assumed to be of the 
form 

(Vu »i, Pi) = %/„„ Xns, P ns ) sin (nAx) [33] 

The coefficients are evaluated in the usual manner. The 
results is a set of coupled ordinary differential equations 
for the amplitude functions and pressure 

dy ns y Q 2 f 2y Q 3 
di t| L 3-n 2 

+ 2 ^mnd/ms ~ 9ms) ~ X n3 

771 = 1 
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with 



r(-ir +n (-ir- n i 

^„ - m + [37] 

Lm+n m-nj 

Because of the rotation terms, all of the equations for the 
coefficients are coupled. However, as is apparent from the 
numerical solution, this coupling is week, and only the 
first few terms must be computed. 

If the capillary number is sufficiently small and the 
mean thickness much less than the wavelength of the to- 
pography (t/oA « 1), the pressure condition, Eq. [36], re- 
duces to 



Pns = 7(TlA) 2 2/n 



[38] 



This is the condition used by Overdiep (10-11). 

If the capillary number were identically zero, Eq. [34] for 
the height fluctuation becomes an algebraic expression 
which can be solved to give 



[39] 



From this result, it is seen that a planar surface (y ns = 0) is 
achieved when the residual solvent fluctuation is zero, 
when the surface tension is very large (7 — > «•), or when the 
surface tension is independent of the composition (r = 0). 
When the pitch is large compared to the thickness 
(2/oA « 1), increasing the mean thickness or decreasing 
the wavelength makes the film more planar. For some con- 
ditions of topography and thickness (VoA > V5), the sur- 
face can "invert": the maximum in the film amplitude cor- 
responds to a minimum in the substrate amplitude. 

A conformal surface (y na = 1) is obtained when there is 
residual solvent remaining and surface tension variations 
with composition are important. If there are no residual 
solvent fluctuations, i.e., if x ns goes to zero, then the surface 
always becomes planar at the steady state. However, in 
many situations of interest, steady state for the fluctua- 
tions occurs on a much longer time scale than for the pri- 
mary coating process. In addition, if the viscosity of the 
material approaches very large values when sufficient sol- 
vent is removed, the fluctuations may become "frozen in" 
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during the spinning. In this case, the product (Ca -n.) is fi- 
nite, and Eq. [38] no longer holds. This appears to be the 
situation for all real films. 

Results 

Parameters. — In order to estimate reasonable values for 
the various parameters needed in the solution to the equa- 
tions, the characteristic length and time scale must first be 
found. In the calculations for the film height on a smooth 
substrate, the final dimensionless thickness is on the order 
of 0.1 (9). Since in most coating operations, the actual film 
thickness is approximately 1 u*m, the characteristic length 
scale must be on the order of 10 u,m. Likewise, the theory 
predicts a dimensionless time of about 3 in order to 
achieve a steady film thickness. The actual time found in 
spinning is about 30s, giving a characteristic time of about 
10s. If the initial fluid viscosity is on the order of 0.01 Pa s 
(10 cp), and the surface tension is 3 x 10~ 4 kg/s (30 dynes/ 
cm), the capillary number is in the range 10** 4 to 10" 5 . For 
surface structures with a period of 1 to 10 p.m, the wave- 
number A is on the order 0.1 to 10. 

Previously (9) two different viscosity models for the so- 
lution were employed. Similar models are used here as 
well 
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' Model 1: t» = {1 + [A(l - x)] n }/{l + [A(l - x°)] n \ [40] 

Model 2: <n = {[1 - p(l - x°)]/[l - Pd - [41] 

Model 1 is a slight variation of the one employed previ- 
ously (9). 

For coating on a smooth substrate, the two models give 
very similar results, despite the fact that the viscosity pre- 
dictions can be orders of magnitude different for small 
values of x. In particular, model 2 predicts infinite viscos- 
ity at a solvent fraction of (1 - 1/0), called the critical sol- 
vent fraction. Presently, there is no reliable model for the 
viscosity over the whole range of concentrations encoun- 
tered. However, Fig. 2 shows some data on standard 
photoresists as a function of solvent content (15). These 
data are clearly consistent with model 2 with a critical sol- 
vent fraction of 0.33. The results are the same for the two 
different manufacturers represented in the figure. This is 
not unreasonable since the solids and solvent are quite 
similar for these positive resists. Despite the better agree- 
ment with model 2, it is instructive to employ both models 
because of the very different types of viscosity predictions. 

Likewise, there are no reliable models for surface ten- 
sion variation with concentration. Photoresists usually 
contain surface active compounds as "leveling agents," 
and so may not show a linear variation. For want of a better 
model, we use the same variation as Eq. [40] 



1.40 



7 = 1 + [B(l - x)] n 7{l + [B(l - x°)] n '} 



[42] 



In addition to the physical parameters needed to solve 
the equations, the initial conditions are also required. Pre- 
viously (9) it was found that the initial film height had no 
effect on the final predicted height, provided the initial 
value was sufficiently large (greater than about 4). For all 
the present calculations, an initial height of 10 was used. 
Numerical experiments indicated that the calculated fluc- 
tuations y i and Xi are also insensitive to their initial values. 
In all the calculations presented here, zero initial fluctua- 
tions are employed. 

Model predictions.— As a base-line case, viscosity model 1 
(Eq. [40]) was used, with a value of A of 1, and an exponent 
of 4. The surface tension model was used with B = 1 and an 
exponent of 2. The surface is assumed to be a simple sinu- 
soid (92s, etc. = 0). The capillary number is 10 " 5 , and A = 1. 
Figure 3 shows the mean height, y ot and the first Fourier 
coefficient of the height fluctuation (the height amplitude), 
y lBy as a function of time with capillary number as a param- 
eter. It is apparent that the two height values can have very 
different time scales. A steady mean height is achieved 
after a dimensionless time of about 2. However, at the time 
the amplitude of the surface wave indicates essentially a 
conformal coating for all values of the capillary number. 
At sufficiently long times, however, the surface is pre- 
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Fig. 6. Effect of surface tension gradient on the surface amplitude, 
n' = 0 (O); 0.1 (O), 0.3 (A), and 1.0 (□). 

dieted to assume a planar value, in agreement with 
Eq. [39]. The time required increases as the capillary num- 
ber increases. For this model, the viscosity of the fluid is al- 
ways finite, so that on a very large time scale, the surface 
tension-driven flow can achieve planarization. The actual 
shape of the surface depends on all of the Fourier coeffi- 
cients, y , a , t/2s, etc. For a sinusoidal substrate, the higher co- 
efficients, y2s, etc., are quite small, indicating an essentially 
sinusoidal surface. 

The effect of the pitch of the topography, A, is illustrated 
in Fig. 4. As the pitch becomes smaller (and hence A 
larger), the amplitude decreases. As observed in the litera- 
ture (3), the finer the structure of the surface, the more 
nearly planar is the coating. Figure 4 also illustrates the 
"inversion" predicted by Eq. [48] for the largest value of A. 
For all but the most widely separated features (A = 0.1), a 
planar coating is predicted at long times. Figure 5 shows 
the effect of residual solvent. For this viscosity model, sol- 
vent remaining in the film has little effect on the type of 
surface which is achieved. The influence of initial solids 
content is also negligible. This is probably more an artifact 
of the viscosity model than a real result. 

The influence of the surface tension gradient is shown in 
Fig. 6. For a constant surface tension (n' = 0 in Eq. [41]), 
the surface is essentially always planar. However, even 
small variations in the surface tension with concentration 




Fig. 5. Effect of residual solvent on the mean height and surface am- Fig. 7. Effect of initial solvent content with viscosity model 2 on 

plitude (y !t ). x R = 0.1 (O; ); x R = 0.3 (A; ). Lines without sym- mean thickness and fluctuations. x 0 = 0.9 ( ; O), 07 (A), and 0.6 

bols refer to mean height. ( — ; A). Lines without symbols refer to mean values. 
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Fig. 8. Effect of residual solvent on the height fluctuations with vis- 
cosity model 2. x R = 0.1 (O), 0.2 (□), and 0.3 (A). If the viscosity re- 
mains finite during the spinning, a planar layer can be obtained. An 
equivalent result is found by changing the critical solvent fraction to 
0.05 (p= 1.05, x R = 0.1; ). 



can have significant effects on the magnitude of the fluc- 
tuations. This is in agreement with Eq. [38] and Overdiep's 
conclusion that high surface tension gradients promote a 
conform al surface. Figure 6 also reveals, as in the previous 
results, that the surface eventually planarizes at a suffi- 
ciently long time. This is due to the finite ultimate value of 
the viscosity. 

As indicated in Fig. 2, viscosity model 2 predicts an un- 
bounded viscosity as the solvent content reaches some 
critical value. Figure 7 shows the predicted mean height 
and surface amplitude for this model. For these calcula- 
tions, a p value of 1.25 is used. Hence, the viscosity be- 
comes unbounded at a solids content of 0.8. The behavior 
of the mean height is very similar to that predicted by 
model 1, as is the value of the height itself. However, if 
there is no residual solvent, the fluctuations become "fro- 
zen in." The viscosity becomes so large that the fluid can- 
not flow, preventing the film from becoming planar. 
Changes in the initial liquid content can modify some- 
what the magnitude of the fluctuations, but a planar sur- 
face can never be achieved. In this case, increasing the ini- 
tial solids content appears to reduce the fluctuations 
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Fig. 10. First Fourier coefficient of the surface fluctuation for cover- 
age of a step function substrate for viscosity model 2 with best fit pa- 
rameters. x 8 = 0.3; x 0 = 07. A = 0.1 (A) # 1.0 (O), and 5.0 (□). 

because the viscosity becomes unbounded before the sur- 
face tension gradients can cause the amplitude to grow 
very large. 

The only way to obtain a planar surface with this model 
is to retain more solvent in the film or to change the critical 
solids content. The predictions are illustrated in Fig. 8. By 
increasing the fraction of retained solvent from 0.1 to 0.3, 
planarization may be achieved using the present parame- 
ters. Similarly, if the critical solids fraction is increased to 
0.95 (a £ value of 1.05), the surface would again be planar. 
The values of x R and p cannot be changed by altering the 
process parameters, however. These depend on the inter- 
actions of the polymer and solvent which make up the re- 
sist and can only be altered by changing the formulation. 

The predicted profiles for a "real" resist system are 
shown in Fig. 9. An initial solvent fraction of 0.7 is used 
with a residual fraction of 0.3. Both viscosity models are 
employed with the "best fit" parameters from the viscosity 
data of Fig. 2. As expected, the mean height for the two 





Fig. 1 1 . The approximation of a square wave substrate with 4 terms 
in the series of Eq. [43] shown in the lower curve. The upper curves rep- 
resent the predicted topography for model 2 viscosity function and best 

fit parameters for A = 0,1 ( ), 1.0 ( ) and 5.0 ( — ). The di- 

mensionless step height, a, is assumed to be 0.1 . 
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models is similar, but the surface fluctuations display 
qualitatively different behavior. Model 1 indicates that the 
surface may be planarized by spinning for a sufficiently 
long time, while model 2 shows that the magnitude of the 
surface amplitude reaches some value and becomes 
frozen. 

Thus far, the substrate has been assumed to be a sinu- 
soid. For the more likely case of a square wave, the to- 
pography function, Eq. [31], is 

9 = sin(Ax) + (1/3) sin (3A X ) + (1/5) sin (5A X ) + . . - [43] 

Calculations were performed using four terms in the 
series. Figure 10 shows the value of y Hf the first term in the 
series for the surface, for three values of the substrate 
pitch and the "real" resist. With features close together 
(large A), the surface tends to be more planar. When the 
features are far apart (small A), the surface is conformal. 
Figure 11 shows the calculated actual substrate and sur- 
face topographies for the same three values of the pitch as 
used in Fig. 10. The first value of the substrate amplitude, 
y l8 , is seen to be a good indicator of the actual nature of the 
surface. 

As mentioned above, the same equations, with some 
slight modifications, may be used to calculate the surface 
topography after a post-apply bake. The results of such a 
computation, for two different spinning times, are shown 
in Fig. 12. In one case, the film is spun for the full time, re- 
sulting in a high viscosity film with a low residual solvent 
content. In the second case, the spin time is much shorter, 
and the film still retains significant solvent. These pre- 
dictions are based on the results for the "real" resist sys- 
tem. As indicated in Eq. [16] the length scale for the two 
problems is the film thickness after spinning. The critical 
solvent content at which infinite viscosity is reached is as- 
sumed to decrease to 0.1 as a result of the high tempera- 
ture. This is the same as the residual solvent fraction. 

The major difference between the two cases is the value 
of the capillary number. For the fully spun film, the viscos- 
ity will be quite high, despite the lowering somewhat by 
the high temperature of the bake. The value of Ca is esti- 
mated to be 0.1. The surface tension is assumed to be con- 
stant for this film. As indicated in Fig. 12, the fluctuations 
change only slightly during the bake. This is a result of the 
large capillary number. For the "wet" film, the viscosity is 
comparatively low. The value of Ca is taken to be 10" 5 . In 
this case, surface tension variations increase the value of 
the fluctuations to some steady value. Decreasing the 



critical solids content to 0.8 reduces the surface amplitude 
somewhat, as shown in Fig. 12. While the theory predicts 
that baking a "wet" film will result in a more planar sur- 
face than when a "dry" film is used, the changes are not as 
dramatic as indicated in Fig. 12. This is because the length 
scale for the two problems is different since the initial film 
thicknesses are different. The actual decrease in the height 
amplitude as a result of partial spinning is about 15%. In 
this case, evaporation of the solvent from the "wet" film 
plays a major role in defining the film topography. 

Conclusions 

A theoretical model for spin coating on topography has 
been developed which relates the shape of the surface to 
the properties of the fluid being spun, as well as to the 
shape of the substrate. Unlike the model for the mean 
thickness of the film, the surface shape is a strong function 
of the physical properties of the material. If the viscosity of 
the fluid can approach a very large value for a finite solvent 
fraction, the shape of the surface becomes frozen during 
the spinning. Planarizing the surface would require de- 
creasing the value of this solvent fraction or increasing the 
amount of residual solvent in the film. In both cases, the 
mobility of the resist would increase. However, changing 
either of these values would require changes to the resist 
formulation. Both of these parameters must be regarded as 
empirical in the present formulation of the problem. The 
value of the critical solvent fraction is related to the chemi- 
cal interaction of the solvent and the solids in the resist. 
The residual solvent fraction is also a result of this interac- 
tion. As such, they are outside the scope of the present 
fluid mechanical study. 

The film surface may be "annealed" somewhat by a post- 
apply bake. The annealing effect is increased by spinning 
only a short time prior to the post-apply bake. In this way, 
more solvent remains in the film at the end of the spin 
cycle allowing more flow during the bake. 

The model predicts that the surface follows the substrate 
for substrates which vary over large distances, and be- 
comes planar for substrates with a short pitch. This is in 
agreement with available experimental data. The model is 
easily extended to more complex substrate geometries by 
expanding the substrate topography in a set of cosine as 
well as sine functions. In that case, an additional set of 
equations for the thickness and composition fluctuations 
must be included in the solution. The calculation for iso- 
lated features would require a different though related ap- 
proach based on Fourier integrals. 

In order to solve the equations, the dimensionless ampli- 
tude of the substrate topography is assumed to be "small." 
As is the case with all perturbation solutions, it is impossi- 
ble to specify how small the neglected terms must be for 
the solution to be valid. If we assume a = 0.1, and use a di- 
mensionless mean height of 0.5, as indicated in Fig. 9, we 
can estimate a ratio of step height to film thickness (cip//t) of 
about 0.2. In many real applications, this ratio may be 0.3 
to 0.5. We would expect the predictions to be qualitatively 
correct for these values, but the actual magnitudes of the 
fluctuations to be different. In particular, the conclusions 
regarding the effects of substrate pitch, retained solvent, 
critical solvent fraction, and post-apply bake should re- 
main unaltered. 

The most serious shortcoming of the model appears to 
be the assumption of uniform properties throughout the 
depth of the film. This assumption requires the inclusion 
of the residual solvent fraction in the model. An improved 
approach would include concentration variations in the 
film, with a concentration dependent diffusivity as well. 

Manuscript submitted June 3, 1988; revised manuscript 
received Jan. 17, 1989. 

LIST OF SYMBOLS 

A parameter in viscosity model 1 

a parameter in evaporation rate 

B parameter in surface tension*model 

^mn constants in equations for the amplitude func- 

tions 

Ca capillary number 
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♦9" differential operator 

e evaporation rate 

h film thickness 

K equilibrium coefficient for the solvent 

Jc mass transfer coefficient 

L characteristic distance 

n exponent in viscosity model 1 

n' exponent in surface tension model 

N variation of the viscosity with composition 

P dimensionless pressure 

p pressure 

<j flow rate 

& rotation parameter 

r radial position 

t time 

V characteristic time 

u dimensionless velocity in z direction 

v velocity 

x solvent fraction 

y dimensionless film thickness 

z direction normal to the substrate surface 

Greek Symbols 

a dimensionless amplitude of surface roughness 

p parameter in viscosity model 2 

T variation of surface tension with composition 

y dimensionless surface tension 

5 exponent in viscosity model 2 

il dimensionless viscosity 

A dimensionless wavenumber of the substrate 

X position 

8 azimuthal position 

M- viscosity 

p density 

£ dimensionless radial position 

a surface tension 

t dimensionless time 

tp substrate topography function 



4 dimensionless position normal to substrate 

ft spin speed 

Subscripts 

B during post-apply bake 

R residual 

s coefficients of the sine terms in the solution 

0, 1, 2, etc. index of coefficients in series solution 

Superscripts 

o initial value 
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ABSTRACT 

Strain distribution in a Si/Si Ju .Ge a strained layer superlattice on Si can be adjusted by epitaxial insertion of a strain re- 
lieved homogeneous Si^Ge^ buffer layer between the substrate and the superlattice. We report on the strain relief in 
Sii. y Ge„ alloy buffer layers with y = 0.225 as a function of layer thickness. The layers are grown by molecular beam epitaxy 
on (100) Si substrates at 450°C. X-ray diffraction, transmission electron microscopy, and Rutherford backscattering are ap- 
plied to the film analysis. Above a critical thickness, i.e., the onset of strain relief, in a first range a strong decrease of film 
strain by formation of misfit dislocations is observed followed by a saturation-like branch with only a weak strain thick- 
ness dependency. A residual strain of about 1/3 of the lattice mismatch is found at a thickness of five times the critical one. 



Presently, strained layer heterostructures and, particu- 
larly, strained layer superlattices (SLSs) are being studied 
for fundamental reasons as well as for device applications. 
Quantum size effects in such artificial, thin layer struc- 
tures combined with the capability of tailoring the elec- 
tronic and optical properties of the individual layers by 
controlling the magnitude of the strain enables a wide vari- 
ety of high performance devices. Recently, the heterosys- 
tem Si-Ge with its relatively large lattice mismatch of 
about 4% has attracted particular attention. Since the early 
report on MBE grown Si/Si^Ge, SLSs on Si substrates by 
Kasper et aL (1) progress in Si MBE technique and devel- 
opment of appropriate doping methods have taken place 
[see for review Ref. (2)]. 



As an interesting result of these developments, Jorke 
et al (3) reported the preparation of n-type modulation 
doped Si/Sio^Ge 05 SLS grown atop an incommensurate 
Sio.75Geo.25 alloy buffer layer on an Si substrate. By using 
the buffer layer both the Si and the Sio. 5 Ge 0 .5 superlattice 
constituents are under strain of opposite sign. This SLS 
showed enhanced room temperature electron mobility by 
formation of a two-dimensional electron gas in the Si lay- 
ers. On the other hand, People et al (4) published a p-type 
modulation doped Si/Si 0 .eGeo. 2 SLS directly grown on an Si 
substrate with a two-dimensional hole gas in the SiGe 
alloy constituents. In this structure only the Sio. B Ge 0 .2 lay- 
ers are strained. The significance of using the buffer layer 
for the band alignment between the Si/Sij^Ge x SLS inter- 



